arXiv: 1507.05848v4 [quant-ph] 6Jun2016 


Generalized Geometric Quantum Speed Limits 


Diego Paiva Pires,'’* Marco Cianciaruso,^’^ '*’^ Lucas C. Celeri,^’* Gerardo Adesso,^’^ and Diogo O. Soares-Pinto'’ “I' 

^Instituto de Fisica de Sdo Carlos, Universidade de Sdo Paulo, CP 369, 13560-970, Sao Carlos, SP, Brasil 
^School of Mathematical Sciences, The University of Nottingham, 

University Park, Nottingham NG7 2RD, United Kingdom 

^Dipartimento di Fisica “E. R. Caianiello”, Universitd degli Studi di Salerno, Via Giovanni Paolo II, 1-84084 Fisciano (FA), Italy 
‘^INFN, Sezione di Napoli, Gruppo Collegato di Salerno, 1-84084 Fisciano (FA), Italy 
^Instituto de Fisica, Universidade Federal de Golds, 74.001-970 Goidnia, Golds, Brazil 

The attempt to gain a theoretical understanding of the concept of time in quantum mechanics has triggered 
significant progress towards the search for faster and more efficient quantum technologies. One of such advances 
consists in the interpretation of the time-energy uncertainty relations as lower bounds for the minimal evolution 
time between two distinguishable states of a quantum system, also known as quantum speed limits. We inves¬ 
tigate how the non uniqueness of a bona fide measure of distinguishability defined on the quantum state space 
affects the quantum speed limits and can be exploited in order to derive improved bounds. Specifically, we estab¬ 
lish an infinite family of quantum speed limits valid for unitary and nonunitary evolutions, based on an elegant 
information geometric formalism. Our work unifies and generalizes existing results on quantum speed limits, 
and provides instances of novel bounds which are tighter than any established one based on the conventional 
quantum Fisher information. We illustrate our findings with relevant examples, demonstrating the importance 
of choosing different information metrics for open system dynamics, as well as clarifying the roles of classical 
populations versus quantum coherences, in the determination and saturation of the speed limits. Our results can 
find applications in the optimization and control of quantum technologies such as quantum computation and 
metrology, and might provide new insights in fundamental investigations of quantum thermodynamics. 

PACS numbers: 03.65.Yz, 03.67.-a 


I. INTRODUCTION 

Quantum mechanics relies on counterintuitive features 
which challenge our merely classical perception of Nature. 
One of the most fundamental quantum aspects lies in the im¬ 
possibility of knowing simultaneously and with certainty two 
incompatible properties of a quantum system [1]. Contrarily 
to the well understood uncertainty relation between any two 
non-commuting observables, the time-energy uncertainty re¬ 
lation still represents a controversial issue [2], although the 
last decades witnessed several attempts towards its explana¬ 
tion [3]. This effort led to the interpretation of the time-energy 
uncertainty relation as a so-called quantum speed limit (QSL), 
i.e. the ultimate bound imposed by quantum mechanics on the 
minimal evolution time between two distinguishable states of 
a system [4—41]. QSLs have been widely investigated within 
the quantum information setting, since their understanding of¬ 
fers a route to design faster and optimized information pro¬ 
cessing devices [42], thus attracting constant interest in quan¬ 
tum optimal control, quantum metrology [43], quantum com¬ 
putation and communication [44]. Interestingly, it has been 
recently recognized that QSLs play a fundamental role also in 
quantum thermodynamics [45]. 
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In a seminal work, Mandelstam and Tamm (MT) [4] re¬ 
ported a QSL for a quantum system that evolves between two 
distinguishable pure states, |i/t(0)) and |i/t(T)), via a unitary 
dynamics generated by a time independent Hamiltonian H. 
The ensuing lower bound on the evolution time is given by 
T > /iarccos(Ki^(T)|i/f(0)>|)/A£', where (AE)^ = {H^} - {H}^ 
is the variance of the energy of the system with respect to the 
initial state. Several years later, Anandan and Aharanov [7] 
extended the MT bound to time dependent Hamiltonians by 
using a geometric approach which exploits the Fubini-Study 
metric defined on the space of quantum pure states. Specif¬ 
ically, they simply used the fact that the geodesic length be¬ 
tween two distinguishable pure states according to the Fubini- 
Study metric, i.e. their Bures angle, is a lower bound to the 
length of any path connecting the same states. Over half a 
century after the MT result, Margolus and Levitin (ML) [16] 
provided a different QSL on the time evolution of a closed 
system whose Hamiltonian is time independent and evolv¬ 
ing between two orthogonal pure states. This bound reads 
as T > jifiKlE), where E - {H) is the mean energy. Al¬ 
though the ML bound is tight, it does not recover the MT one 
whatsoever. Therefore, the quantum speed limit for unitary 
dynamics, when restricting to orthogonal pure states, can be 
made tighter by combining these two independent results as 
T > max(7rIt/(2A£’),;7r^/(2£’)] [17]. 

All these results attracted a considerable interest in the sub¬ 
ject. Giovannetti et al. [18] extended the ML QSL to the 
case of arbitrary mixed states and also showed that entangle¬ 
ment can speed up the dynamical evolution of a closed com¬ 
posite system. A plethora of other extensions and applica¬ 
tions of QSLs for unitary processes has been investigated in 
Refs. [5, 6, 8-15, 19-30]. For example, in Ref. [29] some 
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of us have recently shown that the rate of change of the dis- 
tinguishability between the initial and the evolved state of a 
closed quantum system can provide a lower bound for an in¬ 
dicator of quantum coherence based on the Wigner-Yanase 
information between the evolved state and the Hamiltonian 
generating the evolution. 

Since any information processing device is inevitably sub¬ 
ject to environmental noise, QSLs have been also investigated 
in the nonunitary realm. Taddei et al. [3 1 ] and del Campo et 
al. [32] were the first to extend the MT bound to any phys¬ 
ical process, being it unitary or not. Specifically, Ref. [31] 
exploits the quantum Fisher information metric on the whole 
quantum state space and represents a natural extension of the 
idea used in Ref. [7], whereas Ref. [32] exploits the relative 
purity. Then, Deffner and Lutz [33] extended the ML bound 
to open quantum systems by adopting again a geometric ap¬ 
proach using the Bures angle. These authors have also intro¬ 
duced a new sort of bound, which is tighter than both the ML 
and MT ones, and shown that non-Markovianity can speed up 
the quantum evolution. Some other works have then provided 
a QSL for open system dynamics by using the relative purity, 
whose usefulness ranges from thermalization phenomena [34] 
to the relativistic effects on the QSL [35]. Further develop¬ 
ments include the role of entanglement in QSL for open dy¬ 
namics [36-38], QSL in the one-dimensional perfect quantum 
state transfer [44], and the experimental realizability of mea¬ 
suring QSLs through interferometry devices [39]. Finally, a 
subtle connection was recently highlighted between QSLs and 
the maximum interaction speed in quantum spin systems [40], 
with implications for quantum error correction and the relax¬ 
ation time of many-body systems [41]. 

Distinguishing between two states of a system being de¬ 
scribed by a probabilistic model stands as the paradigmatic 
task of information theory. Information geometry, in par¬ 
ticular, applies methods of differential geometry in order to 
achieve this task [46]. The set of states of both classical and 
quantum systems is indeed a Riemannian manifold, that is the 
set of probability distributions over the system phase space 
and the set of density operators over the system Hilbert space, 
respectively. Therefore it seems natural to use any of the pos¬ 
sible Riemannian metrics defined on such sets of states in or¬ 
der to distinguish any two of its points. However, it is also 
natural to assume that for a Riemannian metric to be bona 
fide in quantifying the distinguishability between two states, 
it must be contractive under the physical maps that represent 
the mathematical counterpart of noise, i.e. stochastic maps in 
the classical settings and completely positive trace preserving 
maps in quantum. Interestingly, Cencov’s theorem states that 
the Fisher information metric is the only Riemannian metric 
on the set of probability distributions that is contractive under 
stochastic maps [47], thus leaving us with only one choice of 
bona fide Riemannian geometric measure of distinguishabil¬ 
ity within the classical setting. On the contrary, it turns out 
that the quantum Fisher information metric [48, 49] is not the 
only contractive Riemannian metric on the set of density op¬ 
erators, but rather there exists an infinite family of such met¬ 
rics [50], as characterized by the Morozova, Cencov and Petz 
theorem [51, 52]. 



FIG. 1. (Color online). Illustration of geometric quantum speed 
limits. The dashed blue curve is the path y in the quantum state 
space representing a generic evolution between an initial state po and 
a final state pj, parameterized by time t e [0, t]. Given a metric 
on the quantum state space, the length of this path is denoted by 
Q(Po,Pt)- The solid red curve denotes the geodesic connecting po to 
Pr, whose length is given by -C(po,Pt)- Quantum speed limits orig¬ 
inate from the fact that the geodesic amounts to the path of shortest 
length among all physical evolutions between the given initial and 
final states: T(po,Pr) ^ Q(Po,Pt) Vy. Such inequality can be inter¬ 
preted as follows. For any given physical evolution y from po to pr, 
and according to any valid metric, the maximum distance between 
the initial po and the final state p^ is the length of the path f^(po,p.r) 
followed by the system. The ensuing minimal time necessary for this 
distance to reach a chosen value is the time at which the path length 
reaches this value. This interpretation provides a neat criterion for 
the saturation of the lower bound on the evolution time, that is when 
the dynamical evolution coincides with a geodesic of the considered 
metric. Here we establish a general family of geometric quantum 
speed limits with respect to an infinite hierarchy of contractive Rie¬ 
mannian metrics on the space of quantum states, unifying and ex¬ 
tending previous results under an information geometry framework. 


In this paper, we construct a new fundamental family of 
geometric QSLs (see Fig. 1) which is in one to one corre¬ 
spondence with the family of contractive Riemannian metrics 
characterized by the Morozova, Cencov and Petz theorem. We 
demonstrate how such non uniqueness of a bona fide measure 
of distinguishability defined on the quantum state space af¬ 
fects the QSLs and can be exploited in order to look for tighter 
bounds. Our approach is versatile enough to provide a unified 
picture, encompassing both unitary and nonunitary dynamics, 
and is easy to handle, requiring solely the spectral decompo¬ 
sition of the evolved state. This family of bounds is naturally 
tailored to the general case of initial mixed states and clearly 
separates the contribution of the populations of the evolved 
state and the coherences of its time variation, thus clarifying 
their individual role in driving the evolution. 

We formulate in rigorous terms the problem of identifying 
the tightest bound within our family for any given dynamics. 
While such a problem is unfeasibly hard to address in gen¬ 
eral, we establish concrete steps towards its solution in prac¬ 
tical scenarios. Specifically, we show explicit instances of 
QSLs which make use of some particular contractive Rieman- 
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nian metric such as the Wigner-Yanase skew information and 
can be provably tighter than the corresponding QSLs obtained 
with the conventional quantum Fisher information. These in¬ 
stances are relevant in metrological settings. Overall this work 
provides one of the most comprehensive and powerful ap¬ 
proaches to QSLs, with potential impact on the characteriza¬ 
tion and control of quantum technologies. 

The paper is organized as follows. In Sec. II we review the 
relation between statistical distinguishability and the contrac¬ 
tive Riemannian metrics on the quantum state space character¬ 
ized by the Morozova, Cencov and Petz theorem. Section III 
provides a new generalized geometric derivation of a family 
of QSLs which is in one to one correspondence with the fam¬ 
ily of such metrics. In Sec. IV we illustrate and compare the 
obtained bounds for both unitary and nonunitary evolutions. 
Finally, in Sec. V we present our conclusions. 


II. GEOMETRIC MEASURES OF DISTINGUISHABILITY 


According to the standard formulation of quantum mechan¬ 
ics, any quantum system is associated with a Hilbert space 
"H and its states are represented by the Riemannian mani¬ 
fold S = Oi'H) of density operators over "K, i.e. the set of 
positive semi-definite and trace one operators over the car¬ 
rier Hilbert space. A Riemannian metric over S is said to 
be contractive if the corresponding geodesic distance X con¬ 
tracts under physical maps, which means satisfy the inequal¬ 
ity X(A(p), A(cr)) < X(p, cr) for any completely positive trace 
preserving map A and any p,cr e S. The Morozova, Cencov 
and Petz theorem provides us with a characterization of such 
metrics in the finite-dimensional case, by constructing a one to 
one correspondence between them and the Morozova-Cencov 
(MC) functions, a function /(f) : R+ ^ R+ which is (/) op¬ 
erator monotone; for any semi-positive definite operators A 
and B such that A < B, then f(A) < f(B)\ (ii) self-inversive: 
it fulfils the functional equation /(f) = tfiljty, and (Hi) nor¬ 
malized: /(I) = 1. Specifically, the Morozova, Cencov and 
Petz theorem states that every contractive Riemannian metric 
assigns, up to a constant factor, the following squared in¬ 
finitesimal length between two neighboring density operators 
p and p + dp [53] 

ds^ gl(dp,dp), (1) 

with 

g^'(A,B)= iTr[Ac^(Lp,Rp)B], (2) 

where A and B are any two traceless hermitian operators, and 


c^(x,y) 


1 

yf(xly) 


(3) 


is a symmetric function, c^(x,y) - c-^(y,x), which fulfills 
c^{ax,ay) = a^^c^(x,y), with /(f) being a MC function, and 
finally Lp,Rp ; Si'H) Si’H) are two linear superopera¬ 
tors defined on the set S('H) of linear operators over "K as 
follows; LpA = pA and RpA = Ap. We stress again that 


each contractive Riemannian metric is arbitrary up to a con¬ 
stant factor. In accordance with Ref. [54], we have chosen the 
factor 1/4 in order to make the entire family of contractive 
Riemannian metrics collapse to the classical Fisher informa¬ 
tion metric when p and dp commute. 

In order to make Eq. (1) more explicit, we can write the den¬ 
sity operator p in its spectral decomposition, p = L 

with 0 < py < 1 and YijPj - ^nd get [54] 


ds^ 


1 

4 



j<i 


(4) 


where dpji {i\dp\l) and we note that the summation is 
constrained to the requirement pj > 0. Equation (4) is cru¬ 
cial since it clearly identifies two separate contributions to 
any contractive Riemannian metric. The first term, which is 
common to all the family, depends only on the populations 
Pj of p and can be seen as the classical Fisher information 
metric at the probability distribution pj. The second term, 
which is responsible for the non uniqueness of a contractive 
Riemmanian metric on the quantum state space, is instead 
only due to the coherences of dp with respect to the eigen- 
basis of p and is a purely quantum contribution expressing the 
non-commutativity between the operators p and p + dp. Fi¬ 
nally, for all the contractive Riemannian metrics that can be 
naturally extended to the boundary of pure states, such that 
/(O) + 0, the Fubini-Study metric appears always to be such 
extension up to a constant factor, so that the non uniqueness 
of a contractive Riemannian geometry can be only witnessed 
when considering quantum mixed states. This is the reason 
for which only the mixed states will be relevant in our anal¬ 
ysis, whose aim is exactly to investigate the freedom in the 
choice of several inequivalent bona fide measures of indistin- 
guishability in order to get tighter QSLs. 

As pointed out by Kubo and Ando [55], among the MC 
functions there exists a minimal one, fmm(t) = 2f/(l -H f), 
and a maximal one, fmaxit) = (1 + 0/2, such that a generic 
MC function f(t) must satisfy fmm(t) < fit) < fmaxit). Inter¬ 
estingly, the maximal MC function is the one corresponding 
to the celebrated quantum Fisher information metric, whereas 
the Wigner-Yanase information metric corresponds to an in¬ 
termediate MC function, fwrit) = (l/4)( V? + 1)^. as illus¬ 
trated in Fig. 2. 

Each of these metrics plays a fundamental role in quantum 
information theory since the corresponding geodesic length 
X, being by construction contractive under quantum stochas¬ 
tic maps, represents a bona fide measure of distinguishability 
over the quantum state space. However, finding such geodesic 
distance is unfortunately a very hard task in general. In fact, 
analytic expressions are known only for the geodesic distance 
related to the quantum Eisher information metric [56], 


X®^(p, cr) = arccos [ -/F(p7o/ ] , (5) 

j is the Uhlmann fidelity, 

and for the one related to the Wigner-Yanase information met¬ 
ric [50], 


where 


Fip,o-) = (Tr ^y[pa-y[p 


cr) - arccos [A(p, cr)] , 


(6) 
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FIG. 2. (Color online). The family of Morozova-Cencov (MC) func¬ 
tions is upper bounded by a maximal one, f„,aAt) = (1 + 0/2, corre¬ 
sponding to the quantum Fisher information metric (red solid line), 
and lower bounded by a minimal one, /„„„(0 = 2f/(l + 1) (purple dot¬ 
ted line). Any MC function f{t) satisfies /m,„(0 ^ /(O ^ fmaxA), i.e., 
its graph falls into the shaded area, as shown in the particular case 
of the MC function fwy{t) = (1/4)(V/+ 1)^ corresponding to the 
Wigner-Yanase information metric (blue dashed line). All the MC 
functions collapse to 1 at r = 1. This reflects the fact that the cor¬ 
responding Riemannian metrics are regular on the maximally mixed 
state. 


where A(p, cr) -TtIsJp yfa) is known as quantum affinity. 


III. GENERALIZED GEOMETRIC QUANTUM SPEED 
LIMITS 


We are now ready to present our main result, that is, a 
family of geometric QSLs which hold for any physical pro¬ 
cess and are in one to one correspondence with the contrac¬ 
tive Riemannian metrics defined on the set of quantum states. 
The most general dynamical evolution of an initial state po 
can be written in the Kraus decomposition as px - £^[po] = 
YjjKjPQK'j\ where {Kj} are operators satisfying Kj - 

I and depending on a set T = {Ai,A 2 ,..., AA of r parameters 
which are encoded into the input state po, in such a way that 
px depends analytically on each parameter (p = 1,..., r). 
In the unitary case, the evolution is given in particular by 
£>a{Pq\ - UxPqU^, where Ux is a multiparameter family of 
unitary operators, fulfilling UxU^ - u\Ux - I. In this case, 
the observables = -ifiUx df,u\, with = dldA/j, are the 
generators driving the dynamics. 

Consider a dynamical evolution px in which the set of pa¬ 
rameters A is changed analytically from the initial configu¬ 
ration A; to the final one Af. Geometrically, this evolution 
draws a path y in the quantum state space connecting px, and 
pXp whose length is given by the line integral - f^ds and 

depends on the chosen metric (see Fig. 1). Since y is an 
arbitrary path between px, and p^^, its length need not be the 
shortest one, which is instead given by the geodesic length 
£j(px,,p.\p) between px, and pxp. Therefore the latter repre¬ 


sents a lower bound for the length of the path drawn by the 
above dynamical evolution. This observation will play a cru¬ 
cial role in the imminent derivation of our family of QSLs, in 
analogy with Refs. [ 7 ] and [31]. 

Since the density operator px evolves analytically with re¬ 
spect to the parameters A, we can write 

r 

dpx = Y^dftpxdA^. (7) 

Let Px - YjjPj\j)(j\ be the spectral decomposition of px, with 
0 < Pj < I and YijPj - 1- We note that both the eigenvalues 
Pj and eigenstates \j} of px may depend on the set of parame¬ 
ters A, i.e. Pi = Pj(A) and \j} = |y(zl)), so that 

d,px = J]^(d,Pj)\j}{j I + Pj[(d,\j)){j I + I)]) , (8) 

j 


and thus 


{j\d^PA\l} = djid^pj + {pi - Pj){j\df,\l), (9) 

where we used the identity = -{jld/ill). Combining 

Eq. (7) and Eq. (9), we get 

r 

{j\dpA\l) = ^ l^jid^iPj + i{pj - pi)^^ji]dAfi , (10) 

where we define s i{j\df^\l). By using Eq. (10), in the case 
of j - I we get 


\{j\dpA\j)f' ^ '^df.pjdyPjdAf.dAy, (11) 

whereas in the case of j + I, by using the fact that dpx is 
hermitian, we obtain 


\{j\dpx\l)\^ = Xi ^PJ - dA,dAy . (12) 

/i,V=I 

Einally, by substituting Eqs. (11) and (12) into Eq. (4), the 
squared infinitesimal length ds^ between px and px + dpx ac¬ 
cording to any contractive Riemannian metric becomes 

r 

ds^ - ^ gfiydAfjdAy , 

/J,V=1 

where 

glv = + Oj^y , 

with 

and 

c^^Pj,Pi){P} - , 

j<i 


(13) 

(14) 

(15) 


(16) 
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referring to, respectively, the contribution of the populations 
of and of the coherences of dpx to the contractive Rieman- 

f 

nian metric tensor gjiy. Herein we restrict to the case where the 
parameters A are time-dependent, (f), for p - 1,..., r, 

and choose the parametrization f e [0, t] —» d (f) such that 
d/ = d(0) and Ap - A(t), where t is the evolution time. 
Now, being the geodesic distance X.Hpo,Pt) between the ini¬ 
tial and final state, po and pj, a lower bound to the length 
({(Pq,Pt) - - fo dt{dsldt) of the path y followed by 

the evolved state p, when going from po to p^, we have 

-C^(Po,Pt) < 40oo,Pt) , (17) 

where 


d-yipO^ Pt) 



f dA^ dAy 
dt dt 


( 18 ) 


Equation (17) represents the anticipated infinite family of 
generalized geometric QSLs and is the central result of this 
paper. Any possible contractive Riemannian metric on the 
quantum state space, and so any possible bona fide geometric 
quantifier of distinguishability between quantum states, gives 
rise to a different QSL. More precisely, we have that both the 
geodesic distance appearing in the left hand side and the quan¬ 
tity ty{p(),pj) being in the right hand side of Eq. (17) depend 
on the chosen contractive Riemannian metric, specified by a 
MC function /. In particular, by restricting to the celebrated 
quantum Eisher information metric, we recover the QSL in¬ 
troduced in Ref. [31]. 

It is intuitively clear that the contractive Riemannian metric 
whose geodesic is most tailored to the given dynamical evo¬ 
lution is the one that gives rise to the tightest lower bound to 
the evolution time as expressed in Eq. (17). In order to deter¬ 
mine how much a certain geometric QSL is saturated, i.e. its 
tightness, we will consider the relative difference 

f 4(PO’Pt)--C^(Po,Pt) 

6y -r-, (19) 

■L'(Pq,Pt) 

that quantihes how much the dynamical evolution y differs 
from a geodesic with respect to the considered metric . 

By minimizing the quantity 6y over all contractive Rieman¬ 
nian metrics, i.e., over all MC functions /, one has a crite¬ 
rion to identify in principle the tightest geometric QSL, of the 
form given in Eq. (17), for any given dynamics y. Eotmally, 
labelling by f* the optimal metric for the dynamics y, the 
tightest possible geometric QSL is therefore defined by 


Nevertheless, in this paper we will move the hrst steps for¬ 
ward towards addressing such general problem, by restricting 
the optimization in Eq. (20) primarily over these two paradig¬ 
matic and physically significant examples of contractive Rie¬ 
mannian metrics, namely the quantum Eisher information and 
the Wigner-Yanase skew information. Quite remarkably, this 
restriction will be enough to reveal how the choice of the 
quantum Eisher information metric, though ubiquitous in the 
existing literature, is only a special case which does not al¬ 
ways provide the tightest lower bound. On the contrary, we 
will show how the Wigner-Yanase skew information metric 
can systematically produce tighter bounds in a number of sit¬ 
uations of practical relevance for quantum information and 
quantum technologies, in particular in open system evolu¬ 
tions. 


IV. EXAMPLES 


In this section we will apply our general formalism to 
present and analyze QSLs based primarily on the quantum 
Eisher information and the Wigner-Yanase skew information 
in a selection of unitary and nonunitary physical processes. 
This will serve the purpose to illustrate how the choice of a 
particular bona fide geometric measure of distinguishability 
on the quantum state space affects the QSLs, therefore provid¬ 
ing a guidance to exploit the freedom in this choice to obtain 
the tightest bounds in practical scenarios. 


A. Unitary dynamics 


We start by restricting ourselves to a closed quantum sys¬ 
tem, so that our initial state po undergoes a unitary evolution 
Pa - UAPqU\. Since the eigenvalues pj of a unitarily evolv¬ 
ing state are constant, dfjpj - 0, we have that 7^^ = 0, and 
thus gfiy - O/^y, along the curve y drawn by the evolved state 
p^. In other words, the coherences of dp a drive the evolu¬ 
tion of a closed quantum system. Moreover, one can easily 
see that - (l/?z){y|A//^|/), where - {Hp, 

(H^) — Tt(paH^), and H'^ - -ihUA dpU\ are the generators of 
the dynamics, so that 


siv 




z 


cHpj,Pi)(Pj - . 


( 21 ) 


{Po,Pt) < iy (po,Pt) , with /* such that inf 6y = 6y , 

( 20 ) 

where the minimization is over all MC functions /. 

Einding this minimum is, however, a formidable problem, 
which is made all the more difficult by the fact that the quan¬ 
tum Eisher information metric and the Wigner-Yanase infor¬ 
mation metric are the only contractive Riemannian metrics 
whose geodesic lengths are analytically known for general dy¬ 
namics (as previously remarked). 


In the following subsections we will focus on, respectively, 
the quantum Eisher information metric and the Wigner- 
Yanase information metric. 


1. Quantum Fisher information metric 

The quantum Eisher information metric corresponds to the 
MC function /(f) = (1 H- f)/2, so that c^(x,y) - 2/(x -i- y) and 





6 


Eq. (21) becomes 


4/( yfx + and Eq. (21) becomes 


^ 2h^ ^ Pj +pi ^ 


( 22 ) 


Moreover, by using the following straightforward inequality 
(Pj - Pif 


Pj + Pi 


Pj + Pi > 


(23) 


we get 


- iZ ^Pj + pi)(MH^,m\^Hi\j} 


jj 


= AHi), 


(24) 


where ^{AH^,AH^) := (l/2)Tr[p.i{A//^, A//(!)] is the sym¬ 
metrized covariance of AH^ and AH^ with respect to the 
evolved state, which reduces to the squared variance of the op¬ 
erator when p - v, i.e. ‘rf(AH^,AH^) - Tr[p,i(A//^)^] = 
By substituting the inequality (24) into 
Eq. (17) we get the new bound 


■CP^(Pq,Pt) 



^{AHf^AHiy 


dA, 


'/i dAy 


dt dt 


(25) 


Although the QSL in Eq. (25) applies to the very general 
r-parameter case, let us restrict for simplicity to the one- 
parameter case where A — t. Consequently, we have that 
^ H, = -ihUidiU] and that the symmetrized covariance 
just reduces to the variance of the observable H, generating 
the dynamics of the system. Therefore, Eq. (25) turns into the 
simpler bound 


WF 

O/iV 



Pj - Pi 
^lp]+ 


2 

I 01A//^^|1)<1|A//;!|;) 


4z^v^- V^^OlA//;]|i)<i|A//;!lj) 

-^Tr([v^,A//;;][Vp,A//,"]) 

YiC^AHf^AH^y), 


(28) 


where C(A//;;,Ai/;!) -(l/2)Tr([^,Ai/;|][ V^,A//;j]) 

reduces to the skew information between the evolved state 
and AH'^ when p - v , C(AH^,AH^) = I(pa,AH^) := 

-(l/2)^T([^/^,AH^f) [57]. By putting Eq. (28) into the 
bound in Eq. (17), we get 


-C'^^(P0,Pr) 


< 



C{AHf„AHi) 


dAj^ dAy 
dt dt 


(29) 


Eor simplicity, let us again analyze the one-parameter case, 
where A - t, ^ H, - -ihUtdiU] and C reduces to 
the skew information I(p,,H,) between the evolved state p, 
and the observable Hr generating the dynamics of the system. 
Therefore, the bound in Eq. (29) turns into 


^Z'^^(Po,Pt) < f dt ^JI(p„H,) 

n tJo 


V2^ 


(30) 


T-^d:^^(pO,Pr) < dt ^{Hf) - <//,)2 

= h-^AE , (26) 

where AE T~^fjdt ^{Hf) - is the mean variance of 
the generator 77,. The following QSL is thus obtained 

T > =2:e^(po,Pr). (27) 

AE 

It is worth emphasizing that the bound in Eq. (27) applies 
to arbitrary initial and hnal mixed states and generic time- 
dependent generators of the dynamics. Moreover, we can im¬ 
mediately see that it exactly coincides with the one reported 
in Ref. [9] and reduces to a MT-like bound, when further re¬ 
stricting to the case of a time-independent generator of the 
dynamics. 

2 . Wigner- Yanase information metric 

The Wigner-Yanase information metric corresponds to the 
MC function /(f) = (l/4)( V? + 1)^, so that {x,y) - 


where we define yflT~^fjdt ^]I{pt. Hr) as the mean skew 
information between the evolved state and the generator of the 
evolution. The QSL thus becomes 

T > — ^!=£^^(po,Pt) ■ (31) 

V2VJ 

As reported by Luo [58], the skew information I(pt, Hr) is up¬ 
per bounded by the variance of the observable 77,, X(p,, 77,) < 

(Hf) - (77,)^, so that VT < AE and 

T > -^£^\p,,Pr) > ^=£.'^^{p,,Pr) . (32) 

V2VJ V2A£ 

The latter QSL strongly resembles the bound expressed in 
Eq. (27) and emerging from the quantum Eisher information 
metric, with the difference lying in the fact that we are now 
adopting the Hellinger angle instead of the Bures angle and 
a V2 factor appears in the denominator. However, when the 
initial and final states commute, we have that the correspond¬ 
ing fidelity and affinity coincide, 7 ’(po,Pt) = A(po,pr), and so 
the Bures angle is equal to the Hellinger angle, which implies 
that in this case the bound emerging from the Wigner-Yanase 
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metric is less tight than the one corresponding to the quantum 
Fisher information by a factor of 1 / V2. 

The above result could be intuitively expected due to the 
strict hierarchy respected by the MC functions corresponding 
to the two adopted metrics. To put such an intuition on rig¬ 
orous grounds, in Appendix A we prove that the geometric 
QSL corresponding to the quantum Fisher information met¬ 
ric, as expressed directly by Eq. (17), is indeed tighter than the 
one corresponding to the Wigner-Yanase information metric, 
when considering any single-qubit unitary dynamics. How¬ 
ever, we leave it as an open question to assess whether this is 
still the case when considering higher dimensional quantum 
systems, or other contractive Riemannian metrics in place of 
the Wigner-Yanase one. 

Quite surprisingly, we will show instead in the next section 
that, for the realistic and more general case of nonunitary dy¬ 
namics, the hierarchy of the MC functions does not automat¬ 
ically translate anymore into a hierarchy of tightness for the 
corresponding QSLs, not even in the case of a single qubit. 
This will reveal original consequences of our analysis in prac¬ 
tically relevant scenarios. 


B. Nonunitary dynamics 

We will now consider two paradigmatic examples of 
nonunitary physical processes acting on a single qubit; de¬ 
phasing and amplitude damping. 


1. Parallel and transversal dephasing channels 

Let us denote by p = (1/2)(I H- r ■ (?) the Bloch sphere 
representation of an arbitrary single qubit state p, where 
r = {r sin cos cf), r sin 6 sin 0, r cos 0} is the Bloch vector, with 
re [0,1], 0 e [0, tt] and 0 6 [0,27r[, while I denotes the 2x2 
identity matrix and <t - {cri, (T 2 , (T 3 ) is the vector of Pauli ma¬ 
trices. 

Let us now consider a noisy evolution of this state governed 
by a master equation of Lindblad form 

^ ^ M’ip) + ^(p), (33) 

where Jff(p) - -i[H,p] describes the unitary evolution gov¬ 
erned by a Hamiltonian H while ^(p) is the Liouvillian 
that describes the noise. We further consider as Hamiltonian 
H - where wq is the unitary frequency, and as Liouvil¬ 

lian 




'-Z 

/=1 


\ 


(34) 


where F is the decoherence rate and o', > 0 with 2/ o'; = 1- 
We can identify two main modalities of dephasing noise. 
When a 3 = 1, the dephasing happens in the same basis as the 
one specifying the Hamiltonian of our system, a case that can 
be referred to as ‘parallel dephasing’. When instead ai - I, 


the dephasing occurs in a basis orthogonal to the one of the 
Hamiltonian, leading to the situation typically referred to as 
‘transversal dephasing’ [59, 60]. We will explore these two 
cases separately. 

a. Parallel dephasing. The parallel dephasing noise lets 
an initial state po evolve as p, - where 

/g-!Wo(/2 Q \ /g-i(Jo(/2 Q \ 

^0 = V^l Q , ?:i = V^l Q _g™of/2l 

(35) 

are the Kraus operators, and q± - (1 + qt)l2 with qt = 
e^^‘ [61]. Notice that the Kraus operators satisfy not only 
Kj - I but also Y}j=t)KjK\ - I, as such a channel is 
unital. The effect of parallel dephasing is exactly the same as 
the one of phase flip and consists in shrinking the Bloch sphere 
onto the z-axis of states diagonal in the computational basis, 
which are instead left invariant. Moreover, cjq describes the 
rotation frequency around the z-axis. One can easily see that 
the evolved state p, has the following spectral decomposition 

Pt ^'^Pj\St,(l)t)j{St,^t\j, (36) 

J=± 

where p± - (1/2)(1 + ro^,) and 
|0r,<^()± = ^[(cos0o±f()|O) + e'<‘"“'+^‘>Vsin0ol^ , (37) 

with = .^cos^Oq + q^sin^Oo and N± a normalization con¬ 
stant. By putting the above equations into Eqs. (15) and (16) 
one obtains, respectively, 

_ rlq^sin^Oo {dq,ldtf 

J — -;;- 


and 



2 2 


cos^Oq {dqtldty 


rgsin^^o c-^(p+, P-). (39) 


The contractive Riemannian metric gf -T + Qf can be in¬ 
terpreted as the speed of evolution of p,. Equation (38), which 
corresponds to the contribution to g^ common to all the MC 
family, is identically zero for all the initial states such that Oq 
is either 0 or n, that are all the incoherent states lying on the 
z-axis of the Bloch sphere (with density matrices diagonal in 
the computational basis), which are indeed left unaffected by 
the parallel dephasing dynamics. Although 7^ is a function 
of the initial purity tq and of time, it does not depend on the 
initial azimuthal angle (pQ since the eigenvalues of the evolved 
state pj do not depend on 0o- Equation (39), which instead 
describes the truly quantum contribution to the speed of evo¬ 
lution gf and depends on the specific choice of the MC func¬ 
tion /, is identically zero for all the incoherent initial states 
such that 00 is either 0 or n. Notice that in the case 0o = 2r/2, 
for initial states lying in the equatorial xy-plane, is nonzero 
only when the frequency uq is also nonzero. Interestingly, 
does not depend on the initial azimuthal angle (pQ as well, even 
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FIG. 3. (Color online). Evolution path lengths {y for parallel dephasing processes by considering the contractive Riemannian metrics cor¬ 
responding to the following MC functions: /gf(0 = fmaxif) (red solid line), fwrit) (blue dashed line) and /„,„(f) (purple dotted line) for (a) 
ro = 1/4, 6o = 7r/4, y8 = 8 (b) ro = 3/4, Og - njA, /? = 8 (c) rg = 1/2, dg = njA, p - 0 and (d) rg = 1/2, dg - njl, = 0. The insets in 
each panel show the relative difference 6y, Eq. (19), for parallel dephasing processes by considering the quantum Eisher information metric 
(red solid line) and the Wigner-Yanase information metric (blue dashed line); such a relative difference can be regarded as an indicator of the 
tightness of the bounds (the smaller Sy, the tighter the bounds). 


though the eigenstates of the evolved state do depend on 00 - 
In summary, the speed of evolution is obviously zero for ini¬ 
tial states belonging to the z-axis, being them invariant under 
parallel dephasing; it is furthermore symmetric with respect 
to the initial azimuthal angle (pQ, and it arises only from the 
populations of the evolving state when starting from the equa¬ 
torial xy-plane with zero frequency wq. 

In Fig. 3 we compare the evolution path lengths appear¬ 
ing in the right hand side of Eq. (17) and corresponding to 
the three paradigmatic examples of contractive Riemannian 
metrics: the quantum Fisher information metric, the Wigner- 
Yanase information metric, and the metric corresponding to 
the minimal MC function. We consider the only initial pa¬ 
rameters that play a role in all the above analysis, i.e. the 
initial purity ro and polar angle Oq, and the dynamical param¬ 
eter /3 = 0 ) 0 /r, while full details on the computation of all the 
quantities appearing in Eq. (17) are deferred to Appendix B. 
First, it can be seen that by fixing the initial purity ro (respec¬ 
tively, polar angle 6o), the speed of evolution increases as we 
increase the initial polar angle 6o (respectively, purity ro). In 
other words, the farther the initial state is from the z-axis (the 


larger is its quantum coherence), the faster the corresponding 
evolution can be. Second, Fig. 3(d) in particular unveils the 
signature of the populations of the evolved state into the speed 
of evolution. Indeed, according to Eq. (39), the purely quan¬ 
tum contribution to the metric is equal to zero for 6o - njl 
and 0)0 — 0 (j3 - 0). Thus, the speed of evolution is de¬ 
scribed solely by the term T given in Eq. (38) and arising 
only from the populations of the evolved state. In this case, 
the speed of evolution remains invariant for any contractive 
Riemannian metric, since T is common to all of them. How¬ 
ever, it is still susceptible to changes depending on the purity 
and time. 

Let us now investigate how the QSLs in Eq. (17) behave 
by considering the quantum Fisher information metric and the 
Wigner-Yanase information metric, whose geodesic lengths 
are known analytically. In the insets of Fig. 3 we compare the 

f 

tightness parameter 6y, as defined in Eq. (19), when consider¬ 
ing these two metrics, for a parallel dephasing dynamical evo¬ 
lution. We can see that for j3 - % the dynamics does not satu¬ 
rate the bound for any of the two metrics, although the quan¬ 
tum Fisher information metric provides in general a slightly 
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FIG. 4. (Color online) Contour plot of the difference Ad-y = 6^'" - 
between the tightness parameter corresponding to the quantum 
Fisher information metric and the one corresponding to the Wigner- 
Yanase information metric, for the parallel dephasing process as a 
function of ro and do, for (a) Ft = 10, wo = 10, and (b) Ft = 10, tuo = 
0.1. The QSL constructed with the Wigner-Yanase skew information 
is tighter than (respectively looser than) the one constructed with the 
quantum Fisher information when ASy > 0 (resp. A6y < 0), as in 
panel b (resp. a). 


tighter QSL. On the other hand, when j3 - 0 and 6q - njl, 
we have that the QSL is saturated for both metrics, whereas 
for = 0 and Oq = njA, it is instead the Wigner-Yanase in¬ 
formation metric that provides us with a slightly tighter lower 
bound. 

More generally, it is sufficient to compare the difference be¬ 
tween the tightness indicators 6^^ - 5^' for the two metrics in 
the whole parameter space of the parallel dephasing model, 
to identify in which regime each of the two corresponding 
bounds is the tightest. This analysis is reported in Fig. 4, 
showing that the Wigner-Yanase information metric does lead 
in general to a tighter QSL when the frequency cjq is suffi¬ 
ciently small. This is in stark contrast with the case of unitary 
evolutions, discussed in the previous section, and constitutes 
a hrst demonstration of the usefulness of our generalized ap¬ 
proach to speed limits in quantum dynamics. 

b. Transversal dephasing. We now focus on the case 
of transversal dephasing noise, which lets an initial state po 
evolve as p, - (Il2)Yj]ijCTipa-j, where S is a 4 x 4 her- 
mitian matrix whose non-vanishing elements are given by 
iS 00 — a + b, S \i — d + f, S 22 — d — f, S 22 — ei — b, Sq 2 ~ 


and S 20 - -ic, with 

fl=i(l+e-“), (40) 

h = cosh (nM/2) , (41) 

2 « e -«/2 

c - —— sinh(QM/2), (42) 

d^^d-e-"), (43) 

/= ^sinh(QM/2), (44) 


where u - Yt, Q. - sjl - and - Wq/F. It is worth¬ 
while noticing that also the transversal dephasing channel is 


unital, i.e. it leaves the maximally mixed state invariant. This 
channel has proven to be of fundamental interest within the 
burgeoning field of noisy quantum metrology, as shown in 
Chaves et al. [59, 60]. More precisely, transversal dephasing 
noise stands as the relevant scenario whereby one can attain 
a precision in the estimation of the parameter cjq that scales 
superclassically with the number of qubits, even if such noise 
applies independently to each qubit (while any superclassical 
advantage is lost in the case of parallel dephasing noise). 

By writing the spectral decomposition of the density oper¬ 
ator p, we get 

Pt , (45) 

where p± - (1/2)(1 -H ro^,) and 

\di, (pt)± = (46) 

^ [[(2a - 1) cos 00 ± |,]|0) + [{b + ic)e‘^° + sin0o|l>] , 
with 

it - -^(2a^M)2cos20o~-i-~^sin^ , (47) 

i, - b^ +c^ + + 2f[b cos(20o) - c sin(20o)] , (48) 

and N± a normalization constant. By putting the above equa¬ 
tions into Eqs. (15) and (16), one obtains expressions which 
are too cumbersome to be reported here. However, when re¬ 
stricting to the relevant case of an initial plus state (which is an 
optimal probe state for frequency estimation), i.e. po = |+)(+| 
with l-H) = (|0) + |1»/V2, one obtains the following simple 
expressions: 

_^ .-W(Qu/2) 

Q4 (i-e-"g)g ’ 

= ^e^"cHpv,P -), (50) 

OB 

where 

& - ^ [cosh (Qu) + Q sinh (Qu) - 4yS^j . (51) 

Let us now analyze the behavior of the QSLs in Eq. (17) 
corresponding to the quantum Eisher information metric and 
the Wigner-Yanase information metric when considering the 
transversal dephasing dynamics. In Fig. 5 we can see that, 
initializing such dynamics with a plus state, it happens that 
for small enough F and oiq the Wigner-Yanase information 
provides a QSL which is tighter (in particular at short times) 
than the one corresponding to the quantum Fisher informa¬ 
tion. One might identify more generally the region of param¬ 
eters in which this behaviour occurs by studying the trade-off 

f 

between the respective tightness indicators 6y for an arbitrary 
initial state, as in the previous case, although such a study does 
not add any further insight and is not reported here. 

Once more, the present analysis shows that our approach 
applies straightforwardly to obtain novel, tighter bounds in 
dynamical cases of interest for quantum technologies, as here 
corroborated in particular for the metrologically relevant case 
of transversal dephasing noise. 
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T 


FIG. 5. (Color online) Plot of the relative difference 6y, Eq. (19), 
indicating the tightness of the QSLs (the smaller 6y, the tighter the 
bounds) corresponding respectively to the quantum Fisher informa¬ 
tion (red solid line) and the Wigner-Yanase skew information (blue 
dashed line), for a qubit initially in the plus state and undergoing a 
transversal dephasing process with (a) F = 0.1 and coo = 0.01; (b) 
F = 0.25 and cug = 0.01; (c) F = 0.25 and wq = 0.033; (d) F = 1 and 

Wo = 1. 


2. Amplitude damping channel 


We now consider another canonical model of noise, namely 
dissipation modelled by an amplitude damping channel acting 
on a single qubit. For the amplitude damping channel we have 
the following Kraus operators 


Ko^ 


0 \ 
VT^j 




(52) 


with A, - 1 - e and 1 /F is the characteristic time of the 
process [61], satisfying only YjjKjKj — I since this chan¬ 
nel is not unital. The effect of amplitude damping consists 
in shrinking the Bloch sphere towards the north pole, or the 
state |0). In this case, it is easy to verify that the evolved state 
pi - TjjKjPoK\ - (1/2)(I + rt ■ a) has the following spectral 
decomposition 


Pt = ^ (I’djidt, 0,1 j , (53) 

f=± 


where p± = 

|(1 + &,) and 


\0t, 4>t)± -- 

1 

[(g,±MO} + e‘•^Oro^|i^ 

/l,sin6»o|l)] , 




(54) 

with 






1 

1 

(55) 



1 - rl -l-/l,(l - rocos0o)^ 

(56) 



A, -H ro(l - /l,)cos0o , 

(57) 


and N± a normalization constant. By putting the above equa¬ 
tions into Eqs. (15) and (16) one obtains, respectively. 




[4', - (1 - T,)(l - rp cos 6>o)^]^ 

16 # 2 ^,( 1 -T,) 


and 


_ r^sin^eo (2 - g,f c^(p+,p-) 
~ 32 # 2(1 - A,) 


(58) 


(59) 


As in the case of the parallel dephasing channel, both contri¬ 
butions !F and to the speed of evolution do not depend 
on the initial azimuthal angle 0o- However, contrarily to the 
parallel dephasing channel case, here the purely quantum con¬ 
tribution vanishes only for 0q = 0, n, whereas the term 'F 
vanishes in neither of these cases nor for Oq = 7r/2, as expected 
due to the fact that now only the north pole, and not the entire 
z-axis of the Bloch sphere, is left invariant by the dynamics. 

In Fig. 6 we compare the evolution path lengths appearing 
in the right hand side of Eq. (17) and corresponding to the 
usual contractive Riemannian metrics, i.e. the quantum Eisher 
information metric, the Wigner-Yanase information metric 
and the metric corresponding to the minimal MC function, by 
changing again the initial purity ro and polar angle 6q. Eirst, 
Eig. 6(e) and Eig. 6(f) exhibit the following behaviour; fixing 
the initial polar angle ffg - 0, the speed of evolution decreases 
as we increase the initial purity tq. This feature highlights the 
fact that the north pole of the Bloch sphere is unaffected by the 
amplitude damping channel. Moreover, according to Eq. (59), 
the purely quantum contribution (3^ vanishes identically for 
00 - 0 and the speed of evolution g and corresponding evolu¬ 
tion path length in Eq. (17) become independent of the choice 
of the MC function /. The nontrivial contribution to the speed 
of evolution is in this case exclusively due to the term 'F which 
depends solely on the populations pj of the evolved state. 

Let us now analyze the behavior of the QSLs in Eq. (17) 
corresponding to the quantum Eisher information metric and 
the Wigner-Yanase information metric (see Appendix C for 
details). In the insets of Eig. 6 we compare the tightness indi¬ 
cators 6y, as defined in Eq. (19), when considering these two 
metrics, for the amplitude damping channel. We can see that 
in this case the Wigner-Yanase information provides a QSL 
which is almost saturated (in particular at short times) whereas 
the quantum Eisher information does not, except in the case of 
00 = 0 where they both realize tight bounds. What is more, in 
Eig. 7 one can see that for almost all initial states, except for 
a small neighbourhood of the north pole (which is the asymp¬ 
totic state of the amplitude damping channel), it happens that 
^QF > ^WY ^ j g j-jjg Wigner-Yanase informa¬ 

tion metric provides us with a dehnitely tighter (and nearly 
saturated) QSL than the quantum Eisher information metric. 

This reveals another important physical mechanism, dis¬ 
tinct from dephasing, in which our generalized analysis leads 
to significantly tighter bounds than those established in previ¬ 
ous literature, in this case clearly demonstrated in almost all 
the parameter space of relevance. This highlights the power 
of our general approach to reach beyond the state of the art. 
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FIG. 6. (Color online). Evolution path lengths iy for amplitude damping processes related to the contractive Riemannian metrics corresponding 
to the following MC functions: /gfit) = fmaxifi (red solid line), /writ) (blue dashed line) and fmm(t) (purple dotted line) for (a) tq = 1/2, 
9o = njl, (b) ro = 3/4, 9q - njl, (c) tq = 1/4, 9a = 7r/4, (d) yq = 1/4, 9a = njl, (e) ro = 1/4, So = 0 and (f) ra = 1/2, 9a = 0. The insets in 
each panel show the relative difference 6y, Eq. (19), for amplitude damping processes by considering the quantum Eisher information metric 
(red solid line) and the Wigner-Yanase information metric (blue dashed line); such a relative difference can be regarded as an indicator of the 
tightness of the bounds (the smaller Sy, the tighter the bounds). 


V. CONCLUSIONS 

Based on the fundamental connection between the geome¬ 
try of quantum states and their statistical distinguishability, we 
have exploited the fact that more than one privileged Rieman¬ 
nian metric appear in quantum mechanics in order to introduce 


a new infinite family of geometric quantum speed limits valid 
for any physical process, being it unitary or not. Specifically, 
each bona fide geometric measure of distinguishability gives 
rise to a different quantum speed limit which is particularly 
tailored to the case of initial mixed states and such that the 
contributions of the populations of the evolved state and of 

































































12 



FIG. 7. (Color online) Contour plot of (a) the tightness indicator 
of the bound specified by the Wigner-Yanase information met¬ 
ric, and (b) the dilference Ady = - 6^^ between the tightness 

parameter corresponding to the quantum Fisher information metric 
and the one corresponding to the Wigner-Yanase information met¬ 
ric, for the amplitude damping process as a function of ro and dg, for 
Ft = 10. The QSL constructed with the Wigner-Yanase skew infor¬ 
mation is nearly globally optimal (as 6^^^ - 0) and tighter than the 
one constructed with the quantum Fisher information (as indicated by 
ASy > 0) in almost the whole parameter space, but for a small region 
in the bottom-right corner (large rg, small Gg, i-S- around the state |0)) 
in which the quantum Fisher information bound is marginally tighter. 


the coherences of its time variation are clearly separated. This 
work provides a comprehensive general framework which in¬ 
corporates previous approaches to quantum speed limits and 
leaves room for novel insights. 

By investigating paradigmatic examples of unitary and 
noisy physical processes and of contractive Riemannian met¬ 
rics, we have seen in fact how the choice of the quantum 
Fisher information, corresponding to an extremal metric and 
being ubiquitous in the existing literature, is only a special 
case which does not always provide the tightest lower bound 
in the realistic case of open system dynamics. In particular, 
for parallel and transversal dephasing, as well as amplitude 
damping dynamics, we defined a tighter quantum speed limit 
by means of another important but significantly less-studied 
Riemannian metric, namely the Wigner-Yanase skew infor¬ 
mation. The bound is useful in practical scenarios of noisy 
quantum metrology, especially in the case of transversal de¬ 
phasing [59, 60]. 

Our unifying approach provides a concrete guidance to se¬ 
lect the most informative metric in order to derive the tight¬ 
est bound for some particular dynamics of interest. We have 
formulated the problem as an optimization of a tightness indi¬ 
cator over all the infinite family of contractive quantum Rie¬ 
mannian metrics. The metric giving rise to the tightest bound 
is identified as the one whose geodesic is most tailored to 
the evolution under consideration, see Eq. (20). While such 
a problem can only be solved in restricted form at present, 
due to the fact that the quantum Fisher information and the 
Wigner-Yanase skew information are the only two metrics ad¬ 
mitting known geodesics, further progress will be achievable 
in case useful advances on the information geometry for other 


relevant metrics are recorded in the future. 

It is important to remark that the family of speed limits pro¬ 
vided in this paper are within the class of MT-like bounds. 
Following [62], it may be possible to implement some ad¬ 
justments to the adopted unified geometric approach in order 
to provide a generalized geometric interpretation to ML-like 
speed limits as well. This will be explored in a further study. 

Our work readily suggests to explore how the non unique¬ 
ness of a contractive Riemannian metric in the quantum state 
space affects also other scenarios of relevance in quantum in¬ 
formation processing. In several of these scenarios, where the 
quantum Fisher information was adopted and privileged, our 
approach could lead to a more general investigation based on 
information geometry. For example, when considering pa¬ 
rameter estimation, one of the paradigmatic tasks of quan¬ 
tum metrology, the inverse of the quantum Fisher informa¬ 
tion metric sets a lower bound to the mean-square error of 
any unbiased estimator for the parameters through the quan¬ 
tum Cramer-Rao bound [49, 63]. This work inspires the quest 
to provide more general bounds on the sensitivity of quantum 
states to evolutions encoding unknown parameters, based on 
the infinite hierarchy of quantum Riemannian contractive ge¬ 
ometries. It is useful to recall here that the Fisher information- 
based quantum Cramer-Rao bound for single parameter esti¬ 
mation can only be achieved asymptotically in the limit of a 
large number of probes, and upon performing an optimal mea¬ 
surement given by projection into the eigenbasis of the sym¬ 
metric logarithmic derivative, which is typically hard to im¬ 
plement in the experimental practice [63]. In the realistic case 
of a finite number of probes, corrections to the bound provide 
tighter estimates to the attainable estimation precision; these 
corrections have been first investigated in Ref. [64] for the 
case of the quantum Fisher information. Motivated by more 
recent works by Brody [65] and [66], in which the Wigner- 
Yanase skew information has been interpreted rather naturally 
as the speed of mixed quantum state evolution, and by the 
analysis of the present work (which includes metrologically 
relevant settings such as frequency estimation under transver¬ 
sal dephasing [59, 60]), we believe it is a worthy outlook to 
investigate finite-size corrections to the Cramer-Rao inequal¬ 
ity based on the Wigner-Yanase information, in order to de¬ 
termine how tight the bound can be for practical purposes, in 
particular for the estimation of parameters encoded in open 
system dynamics. 

Furthermore, within the burgeoning field of quantum ther¬ 
modynamics, our approach could provide an infinite class of 
generalizations of the classical thermodynamical length [67], 
originally based on the unique classical contractive Rieman¬ 
nian metric, to the quantum setting. Again in the context of 
quantum thermodynamics, due to the close connections be¬ 
tween geometry and entropy, it might be interesting to in¬ 
vestigate the role played by the non uniqueness of a contrac¬ 
tive Riemannian geometry on the quantum state space in the 
existence of many second laws of thermodynamics [68]. In 
the study of quantum criticality, within the condensed matter 
realm, a geometric approach based on the fidelity, i.e. on the 
quantum Fisher information metric, proved to be fruitful [69]. 
Along the lines of this work, one could apply more general 
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tools associated with any quantum Riemanniann contractive 
metric, in order to seek for further insights and sharper identi¬ 
fication of quantum critical phenomena. 

Finally, the general approach presented in this paper to pin¬ 
point the tightest speed limits in quantum evolutions is readily 
useful for applications to quantum engineering and quantum 
control. Specifically, the present study allows one to certify 
that, in a particular implementation, quantum states have been 
driven at the ultimate speed limit [42] and their evolution can¬ 
not be sped up further: this occurs whenever saturation of one 
of our bounds is demonstrated. As our various examples show, 
this is not possibly verifiable only by considering the standard 
bound based on quantum Fisher information. For single-qubit 
evolutions, we showed that the latter is in fact the tightest for 
the idealized case of unitary dynamics, while our novel bound 
based on the Wigner-Yanase skew information can instead be 
signihcantly tighter in the most common instances of open dy¬ 
namics, yielding effectively the optimal bound (even among 
all the other unverifiable Riemannian metrics) for amplitude 
damping dynamics, as certified by a nearly vanishing tight¬ 
ness indicator in such case. Given that the Wigner-Yanase 
skew information is experimentally accessible [70], one can 
readily apply our results to current and future demonstrations 
to benchmark optimality of controlled quantum dynamics in 
the presence of such ubiquitous noise mechanisms. 

In this respect, we would like to point out that an exper¬ 
imental investigation of the main results presented here, for 
both closed and open system dynamics, can be achieved in 
particular using a highly controllable Nuclear Magnetic Reso¬ 
nance setup, with no need for a complete quantum state to¬ 
mography. In fact, dephasing and amplitude damping are 
naturally occurring sources of decoherence in such an imple¬ 
mentation, and our results can be accessed by means of spin 
ensemble measurements, which constitute the conventional 
types of detection in such a technique [71, 72]. An experi¬ 
mental investigation as described deserves a study on its own 
and will be reported elsewhere. 


than the one corresponding to the Wigner-Yanase information 
metric, when considering any single-qubit unitary dynamics. 

Let us consider a one-qubit state po - (1/2)(I -t- ro ■ a), 
where ro = {^o sin Oq cos (^q, rg sin ffg sin fo cos ffg) is 

the vector of Pauli matrices, which undergoes the generic uni¬ 
tary evolution y specihed by p, - U,poUj. We want to prove 
that for any po and Ut, the following holds (we will drop the 
subscript y in the remainder of this section for simplicity) 

6^^ < , (Al) 

where 6^ is the tightness indicator corresponding to the con¬ 
tractive Riemannian metric with MC function /, as defined 
in Eq. (19). In order to prove the above inequality, we just 
need to prove that 

{^^(P0,Pr) £^^(pO,Pr) 

('^HpO,Pr)- -C^^(pO,Pr)^ ^ ^ 

where we denote by £^(po,pr) and X.^(po,Pt), respectively, the 
path length along the given unitary evolution and the geodesic 
length between initial and final state, po and pj - UrPoUl, 
according to the contractive Riemannian metric g^. 

We know that 

£^(po,Pr) = £ (A3) 

where -'T + Q^ and 



Qt = ]^cf{p\,P2){p\ - Pif^n^u , (A5) 

with Pj being the eigenvalues of the evolved state p,, J[ji - 
i being quantities that depend on the eigenstates \ j) of 

the evolved state p,, and finally 
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c^'^iPuPi) = —-— = 2 , (A6) 

Pi +P2 

c'^^(PuP2) = - - - -^ ■ (A7) 

(V/h + yfpy 

Since the eigenvalues pi 2 - (l+ro)/2 of the unitarily evolving 
one-qubit state are time independent, it immediately follows 
that 


T = 0, 


c'^^(Pl,P2) = 


1 H- 




(A8) 

(A9) 


Appendix A: Unitary dynamics 


so that 


In this Appendix we prove that the geometric QSL corre- 

^^^(P0,Pr) 

jcS^(pi,P2) A 

1 

+ 

spending to the quantum Eisher information metric is tighter 


SC'^^(P1,P2)~ ' 

2 


(AlO) 
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On the other hand, due to Eqs. (5) and (6), we have that 

£P’"{pQ,Pr) _ arccos[VE(po,PT)] 
arccos[A(po,pT)] 


where the analytical formulae of the fidelity F(pQ,pj) and the 
affinity A(po,Pt) for any pair of one-qubit states po and pr are 
the following [29, 73, 74] 


F(po,Pr) = ^ [l + ih ■ iv + V(l-|^oF)(l-|ivP)](A 12 ) 


A(po,Pr) = 7 


^0 ^0 


-fo-rr 


kol kr 


(A13) 


with 


= Vl + lol ± Vl - \ra\ , (A14) 


and a - 0,T. 

Being p^ = UrPoUl, it happens that 

lihl - livi - ro , (A15) 

ih ■ iv = |ihl|ivl cos((/?o,t) = cos(95o,t) , (A16) 


with ^o,T being the angle between the vectors ro and iy, so that 
we simply get 


£^^(pO,Pr) 

jy^(pO,Pr) 


arccos 


arccos -^5(2 - H- cos(i,c)o,t)) 

{5 cos(,,oo,r)(l - ^1 -^o) + 1 + 


(A17) 


Overall, by collecting Eqs. (AlO) and (A17), we need to 
prove that 


A 


1 + 

2 


< 


(A18) 


arccos -yj 

IT 

{2- rl + rl cos(^o,t)) 


arccos 

If 

cos((,oo,t) I 

1-v 




The difference between the right hand side and left hand 
side of the above inequality is represented in Eig. 8 as a func¬ 
tion of ro 6 [0,1[ and ipo^ e [0,7r]. As it can be easily seen, 
this difference is always non-negative, i.e. the inequality (A2) 
is always satisfied. In particular, this difference is zero when 
either ro = 0 or ro ■ iv = r^, i.e. when ip^ j - 0, as it can be 
proved by checking that 


arccos 


lim 


V 7^^ “ '■q + ^0 ■ A) 


arccos (i (l - + 1 + 


A 


1 H- 




.(A19) 


0 0.05 0.10 0.15 0.20 0.25 



FIG. 8. The difference between the right hand side and left hand side 
of Eq. (A18) as a function of ro and of the relative angle ipo^r- 


Some remarks are now in order. One the one hand, 
Eq. (AlO) can be trivially generalized to any pair of contrac¬ 
tive Riemannian metrics and g* as follows 


{Hpq,Pt) _ \ cHp\,P2) 

{''ipa,PT) ycHpuPi) ' 


(A20) 


However, we also note that this is only true in the one-qubit 
unitary case. If the dimensionality of the system is higher than 
2, then becomes a non trivial sum, as defined in Eq. (16), so 
that the various time-independent coefficients c^ipj,pi) can¬ 
not be extracted from the path integral as we have just done 
above. On the other hand, Eq. (A17) seems hard to generalise 
to other contractive Riemannian metrics, since the analytical 
expressions of the corresponding geodesic lengths are still un¬ 
known. We can thus leave the following conjecture: for any 
one-qubit unitary dynamics, the quantum Eisher information 
is the metric which provides the tightest QSL among all con¬ 
tractive Riemannian metrics on the quantum state space, i.e., 
it is the metric solving the optimization problem in Eq. (20). 
Here we have shown that the conjecture holds when the opti¬ 
mization is restricted to the quantum Eisher information and 
the Wigner-Yanase skew information metrics only. 


Appendix B: Parallel and transversal dephasing 

In this Appendix we show how to compute all the quantities 
playing a role in Eq. (17), for the parallel and transversal de¬ 
phasing dynamics, when considering the quantum Eisher in¬ 
formation metric, Wigner-Yanase information metric, and the 
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metric corresponding to the minimal MC function. 

Let us start from the case of parallel dephasing, and from 
the quantum Fisher information metric. Recall that its MC 
function is the maximal one, i.e. fQpit) — fmax(t) = (1 + 0/2, 
implying that c^(x,y) - 2f(x+y) and, since Tr(p,) = P++P- - 
1, cHp+,P ) - 2. It is possible to provide an analytical ex¬ 
pression for the path length corresponding to such MC func¬ 
tion as follows 

ty^iPo^Pr) = ^ a/^ 1 +^^)<AZ) [£'(arcsinQ',>f2) 

-£(arcsin, (Bl) 

where E(y,K^) - fj dy' -^1 - sin^ y' is the elliptic integral 
of second kind, p - cuolT, + 0^) and 


<AZ) = 1 - r/cos20o , X 1- , (B2) 

\ <AZ) 

with (AZ) = Tr(ptcr^) - Tr(p,crp^ being the variance of the 
Pauli matrix ctj.. Specifically, when considering the case in 
which a>o - 0, i.e. p - 0 and /c = 0, the path length in Eq. (Bl) 
becomes 

ef{po,Pr) = ^ [ arcsin O'- arcsin (ere . (B3) 

Furthermore, the geodesic length corresponding to such 
metric can be readily obtained by using Eqs. (5) and (A12), 
with 

|ivl = ro^T, (B4) 

ro ■ iv = ^0 [cos^0o + qr^m^Oo cos(woT)j . 

Considering now the Wigner-Yanase information metric, 
we know that the corresponding MC function is fwvit) - 
(l/4)( Vf + 1)^ so that cf{p+,pS) - 4/( yfp^ + yfpZf-. Hence 
the path length corresponding to the Wigner-Yanase informa¬ 
tion metric is given by 


c^{p+,p-) - l/(2p+p_) and so its path length has the fol¬ 
lowing analytical form 

^™”(po,Pr) = ^ (arcsina- arcsin . (B8) 

Unfortunately, the geodesic length of the metric correspond¬ 
ing to the minimal MC function is still analytically unknown. 

On the other hand, regarding transversal dephasing noise, 
we are able to compute analytically both the geodesic 
distances corresponding to the quantum Fisher informa¬ 
tion and the Wigner-Yanase information metric by using 
Eqs. (5), (6), (A12) and (A13) with 

livi = roir (B9) 

ro ■ iv = rllila - l)cos^0o + [b + /cos(20o)]sin^0o] ■ 

Notice that we used the constraint a + d - 1, with a and d 
given respectively in Eqs. (40) and (43). The expressions of 
the path lengths are instead too cumbersome to be reported 
here in the transversal dephasing case. 

Appendix C: Amplitude damping 

In this Appendix we show how to compute all the quan¬ 
tities playing a role in Eq. (17), for the amplitude damping 
dynamics, when considering the aforementioned three met¬ 
rics. Similarly to the previous case, we are able to compute 
analytically both the geodesic distances corresponding to the 
quantum Fisher information and the Wigner-Yanase informa¬ 
tion metric by using Eqs. (5), (6), (A12) and (A13) with 

iV'F) = ro {[/Ir + (1 - Ar)ro cos 6»o] cos 6»o -H 01 - Ar rosin^0o} , 

(Cl) 

whereas the geodesic length relative to the minimal MC func¬ 
tion is not known analytically. 

Now we focus on the evaluation of the path length {y{p{),pr) 
for these metrics. For the quantum Fisher information metric, 
the analytical expression of the path length is given by 


£^^(po,Pr) = ^ ^(l +P^){AZ} du 'I'„(ro, 0o, k) 


where 


'Vu{ro,0Q,K) ^ y D.u{rQ,0Q,K) 




1 — K^a^e 


1 - a^e 


2^-2u 


(B5) 


(B6) 


and 


^(po,Pr) = Vl - (e (arcsin m, 

-E (arcsin , 


where 


rhi0eQ 


' 2(1 - rocosffo) 


11 - ro cos 00 
2(1 - £2) 


(C2) 


(C3) 


Q„(ro,6»o,x) 


r2(cos20o + e ^“/<:^si00o) 

1 -I- ^(1 - rQCOs20o)(l - a2e-2ii) 


(B7) 


Furthermore, the geodesic length corresponding to such 
metric can be readily obtained by using Eqs. (6), (A13) 
and (B4). 

Finally, when considering the metric corresponding to the 
minimal MC function, f„m(t) - 2tl(\ + t), we have that 


For the Wigner-Yanase information metric, the path length 
is evaluated as 

f^^iPChPr) = ^ Vl -£2 du ?ire^“/2 g) ^ ((-4) 


where 


^(ro, 6>o, £) 


Q„(ro,6>o,e) 


1 — e^'uPe “ 

1 — tD'2e““ 


(C5) 
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and 


where 


s^[l + 2(1 - e2)OTV“]2 

1^1 +2(1 - e^)TiTe“"/2 Vl — Tjp-e^" j 


Q„(ro,0o,e) = 


■ (C6) 


Finally, for the metric corresponding to the minimal MC 
function, the path length is given by 


0(x) = arctan 


+ ^n 

2 


Va^ + sin^xy 


+ A arctanh (cos x)+ 


2 cos X 


2cos2(x/2) + A (a + Va^ + sin^x j 


(C8) 


^"(Po.Pt) = 0(<l>r) - 0(<I>o) , (C7) O,, = arcsin (nTe-r«/2) and A^ = sV[4(l - e^)]. 
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